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Abstract: The latest CERES FM-5 instrument launched onboard the S-NPP spacecraft will 
use the VIIRS visible radiances from the NASA Eand Product Evaluation and Analysis 
Tool Elements (PEATE) product for retrieving the cloud properties associated with its 
TOA flux measurement. In order for CERES to provide climate quality TOA flux datasets, 
the retrieved cloud properties must be consistent throughout the record, which is dependent 
on the calibration stability of the VIIRS imager. This paper assesses the NASA calibration 
stability of the VIIRS reflective solar bands using the Eibya-4 desert and deep convective 
clouds (DCC). The invariant targets are first evaluated for temporal natural variability. It is 
found for visible (VIS) bands that DCC targets have half of the variability of Eibya-4. 
For the shortwave infrared (SWIR) bands, the desert has less variability. The brief 
VIIRS record and target variability inhibits high confidence in identifying any trends 
that are less than ±0.6%/yr for most VIS bands, and ±2.5%/yr for SWIR bands. None of 
the observed invariant target reflective solar band trends exceeded these trend 
thresholds. Initial assessment results show that the VIIRS data have been consistently 
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calibrated and that the VIIRS instrument stability is similar to or better than the 
MODIS instrument. 

Keywords: satellite calibration; S-NPP VIIRS; radiometric stability; MODIS; CERES; 
invariant calibration targets 


1. Introduction 

The Clouds and the Earth’s Radiant Energy System (CERES) [1] provides the climate community 
with a continuous global top-of-atmosphere (TOA) shortwave and longwave flux observation record 
beginning in 2000 using CERES instruments mounted onboard the Terra, Aqua and Suomi National 
Polar-orbiting Partnership (S-NPP) spacecrafts. CERES relies on the accompanying onboard imager to 
provide the necessary cloud retrievals in order to accurately convert the observed CERES instrument 
footprint radiances into fluxes. The conversion uses Angular Directional Models (ADM), which are 
dependent on solar, viewing and azimuthal angles, surface type, and cloud properties [2,3]. The recorded 
CERES dataset must have consistent cloud properties as not to bias the radiance-to-flux conversion. Any 
drift in the imager calibration may cause spurious cloud retrievals. CERES does not incorporate the 
official NASA Goddard Space Flight Center (GSFC) imager cloud product properties [4], and thus 
applies its own retrieval algorithm, specifically designed to be stable throughout the record, and 
consistent in providing cloud properties for all identified cloud pixels [5]. This method ensures that the 
cloud properties do not introduce any artificial trends in the radiance-to-flux conversion. 

The latest of the CERES instruments. Flight Model-5 (FM-5), was launched aboard the S-NPP 
spacecraft on 28 October 2011 from the Vandenberg Air Force Base, California [6]. Another key 
instrument onboard S-NPP is the Visible Infrared Imaging Radiometer Suite (VIIRS), which collects 
the Earth view data in 22 spectral bands with wavelengths ranging from 0.41 to 12 pm [7,8]. The cloud 
retrievals provided in the CERES FM5 NPP Single Scanner Footprint (SSF) Edition 1 product will be 
based on the simultaneous visible radiance measurements from VIIRS. In order for CERES to provide 
temporally consistent cloud property retrievals, the radiometric stability of the VIIRS visible sensors 
must be assessed and ensured. This study provides an independent assessment of the radiometric 
stability of the VIIRS reflective solar bands calibration using two pseudo-invariant Earth targets: 
Eibya-4 desert in Saharan Africa and deep-convective cloud (DCC) pixels acquired over the Tropical 
Western Pacific (TWP) region. 

2. Sensor Overview and Dataset Selection 

VIIRS is nominally operated in a sun-synchronous orbit at 13:30 local equatorial crossing time. It is 
a scanning radiometer with a large swath width of 3040 km. The ground sampling distance at nadir is 
375 m for high-resolution imagery bands (I-bands) and 750 m for moderate resolution bands 
(M-bands). Because of the better radiometric stability and accuracy, the M-bands are recommended for 
quantitative studies, while the I-bands are useful for producing higher resolution images [8]. Out of the 
22 spectral bands in VIIRS, 14 lie in the visible (VIS), near infrared (NIR), and shortwave infrared 
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(SWIR) regions, and are referred to as the refleetive solar bands (RSB). Seven VIIRS bands are in the 
mid-wave infrared (MWIR) and long-wave infrared (LWIR) regions and, therefore, referred to as the 
thermal emissive bands (TEB). The VIIRS sensor also earries a panehromatie day-night band (DNB). 
The VIIRS RSB are ealibrated on-orbit using a solar diffuser (SD) and a solar diffuser stability 
monitor (SDSM) [7,8]. The SD ealibration data are eolleeted every orbit. The SDSM, eurrently 
operated on a daily basis, uses simultaneous observations of the sun through an attenuator and the 
diffuser to traek SD on-orbit degradation. The ealibration eoeffieients derived from the SD and SDSM 
observations are used to update the VIIRS ealibration lookup tables (LUTs), throughout the mission. 
The VIIRS ealibration LUTs are eurrently updated on monthly basis. 

The CERES projeet does not use the VIIRS level IB radianee sensor data reeord (SDR) produet 
available at the NOAA Comprehensive Large Array-data Stewardship System (CLASS) for eloud 
retrievals, but rather uses an espeeially reformatted subset of the NPP VIIRS AS3I00 dataset [9] 
generated by the NASA Land Produet Evaluation and Analysis Tool Elements (PLATE). The Land 
PLATE uses mission-long, eonsistently reproeessed ealibration LUTs provided by the NASA VIIRS 
Charaeterization Support Team (VCST). This effort enables an independent evaluation of the VIIRS 
SDR quality and supports the seienee diseipline teams in their Environmental Data Reeords (EDRs) 
assessment and validation effort. The NASA Land PLATE ASS 100 subsetted data for the CERES 
projeet are arehived at the Atmospherie Seienee Data Center Distributed Aetive Arehive Center 
(ASDC DAAC), and eontain radianee data for 8 VIIRS RSB. This paper foeuses on the stability 
assessment of these bands only. The eharaeteristies of these bands are listed in Table 1. In this study 
the NIR band is elassified with the VIS bands for eonvenienee. Bands with wavelengths less than 1 pm 
are referred to as VIS bands, and those greater than I pm as SWIR bands. 

Table 1. VIIRS RSB available in the VIIRS Land PLATE ASS 100 subsetted data arehived 

at the ASDC. 


Band 

Central 

Wavelength (nm) 

Spatial Resolution at Nadir (m) 

Speetral Region 

MS 

488 

750 

VIS 

M4 

555 

750 

VIS 

M5 

672 

750 

VIS 

M7 

865 

750 

NIR 

MIO 

1610 

750 

SWIR 

Mil 

2250 

750 

SWIR 

11 

640 

375 

VIS 

IS 

1610 

375 

SWIR 


3. Methodology 

Pseudo-invariant Earth targets have been used extensively to monitor the post-launeh radiometrie 
stability of satellite sensors for over two deeades [10-20]. These targets primarily inelude subtropieal 
deserts, polar iee sheets, and deep eonveetive elouds. The premise of this method is that the TOA 
refleetanee measured over desert and polar-iee seenes, or as in the ease of DCC, the mode and mean 
statisties of eombined DCC pixel response, are effeetively eonstant over time. Any temporal ehanges 
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observed in the sensor response are then eonsidered as instrument degradation. Hereafter, 
pseudo-invariant Earth targets will simply be ealled invariant targets, with the understanding that for 
the purpose of viearious ealibration of satellite sensors, these seenes are eonsidered to be 
radiometrieally and temporally stable. 

3.1. Desert Approach 

Several past researehers have reported on the qualifieation of invariant desert sites for 
ealibration [13,21-23]. Based on these eriteria, the most invariant desert targets are found in the 
Saharan and Arabian Desert regions. The Libya-4 desert (28.6°N latitude and 23.4°E longitude) 
loeated in the Saharan Desert is the most well-known of the six referenee standard test sites 
reeommended by the Committee on Earth Observation Satellites (CEOS) Infrared and Visible Optieal 
Sensors (IVOS) subgroup for monitoring the post-launeh degradation of the spaee-bome optieal 
sensors [24]. Owing to its high surfaee refleetanee, high spatial homogeneity, exeellent radiometrie 
stability (less than 2% in VIS and NIR [21]), low aerosol loading, and minimal eloud eover, Libya-4 is 
eonsidered to be one of the best invariant targets on Earth for post-launeh satellite ealibration and 
validation purpose [20,21,24]. Libya-4 has been used for monitoring the temporal stability of many 
satellite sensors ineluding Landsat Multispeetral Seanner (MSS) [19], Thematie Mapper (TM) [16], 
and Enhaneed Thematie Mapper plus (ETM-I-) [17], AVHRR [12], and Aqua and Terra MODIS [17,18]. 
This study also uses Libya-4 as a radiometrieally invariant target to assess the ealibration stability of 
the 8 VIIRS RSB listed in Table 1. 

The observed pristine elear-sky TOA refleetanee over an invariant desert site is dependent on the 
anisotropie surfaee, the gaseous absorption of water vapor, ozone, and other traee gases, as well as 
aerosol effeets. The measured refleetanee will also eontain the natural variability of the desert surfaee 
refleetanee and any ealibration artifaets of the observing instrument. Sinee anisotropie effeets are 
speetrally dependent, it is preferable to eonstruet a sensor speeifie empirieally derived bidireetional 
refleetanee distribution funetion (BRDF) model [18,20]. Generally, the BRDF is eonstrueted 
using the first few years of data, when the instrument is most understood, and all known ealibration 
anomalies have been removed. The eurrent VIIRS data reeord is too short to derive a robust 
BRDF model. 

For this study, MODIS bands, that are similar to VIIRS bands in Table 1, are used to generate the 
band-speeifie BRDF models. The MODIS Colleetion 6 (C6) radianees offer signifieant improvements 
over Colleetion 5 (C5) in terms of radiometrie stability and aeeuraey [25-27]. At the time of this study 
only the Aqua-MODIS level IB C5 files were available and the C5 radianees were eonverted to C6 
using ealibration eoeffieients provided by the MODIS ealibration support team (MOST). Figure 1 
displays the speetral response funetions (SRF) of the MODIS and VIIRS band-speeifie pairings. The 
0.55 pm, 0.65 pm (with II), and 0.86 pm MODIS and VIIRS bands are very similar, whereas the 
0.48 pm, 0.65 pm (with M5), 2.2 pm bands have little overlap. Aqua is the preferred MODIS 
instrument beeause both Aqua and NPP are in the same 13:30 sun-synehronous orbit, thereby avoiding 
any diurnal BRDF effeets. More than half of the Aqua-MODIS band 6 (1.6 pm) deteetors, however, 
are inoperable [28]. For the VIIRS MIO and 13 bands, the BRDF is, therefore, derived using the 
Terra-MODIS band 6 data. 
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A 0.5° X 0.5° latitude by longitude area centered at 28.6°N latitude and 23.4° E longitude defines 
the Libya-4 region-of-interest (ROI) used in this study. To ensure that only pristine clear-sky 
reflectances are used to assess the stability, a spatial pixel-level homogeneity threshold, as described 
in [20], is used. The standard deviation of the pixel-level reflectances for pristine clear-sky 
observations are consistent and the lowest of the all observations. This threshold is effective in 
removing thin clouds and aerosol events caused by windy conditions. Only near-nadir viewing zenith 
angle (VZA) conditions (VZA < 10°) are considered to minimize the BRDF effects. 


Figure 1. Comparison of the spectral response functions (SRFs) of the matching spectral 
bands of the MODIS (in red) and VIIRS (in green) instruments. Band numbers are posted 
in the upper left corner. 
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Bhatt et al. [29] found that under near-nadir viewing condition, the MODIS visible radiances over 
Libya-4 were stratified by scattering direction or relative azimuth angle (RAA), and can be simply 
modeled with two linear functions of cosine of solar zenith angle (SZA). We used the angular 
convention of CERES ADM [2] that defines the RAA as the azimuth angle position of the satellite 
relative to the solar plane. Using this convention, all measurements with RAA < 90° are considered 
backward scattering, whereas measurements with RAA > 90° are considered forward scattering. 
Figure 2a shows the Aqua-MODIS Band 1 (0.65 pm) TOA radiances plotted as a function of cosine of 
SZA and scattering direction. It clearly shows that the observed TOA radiances are slightly brighter in 
the backward scattering than in the forward scattering direction. The modeled radiances have been 
normalized to an Earth-Sun distance of 1 AU. The MODIS-based forward and backward regression 
models are used to compute the predicted VIIRS TOA radiance for a given SZA and scattering 
direction. The ratio of the observed and predicted VIIRS radiances should be mostly free from angular 
dependencies. This ratio can thus be trended over time to monitor the stability of the VIIRS sensor. 

Figure 2b shows the temporal trend of the VIIRS M5 channel ratio, which exhibits a bias of 5% 
(Aqua-MODIS is darker). A constant bias is observed in the ratio, which is primarily caused by the 
spectral differences between the VIIRS M5 and MODIS Band I as shown in Figure I. The incoming 
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spectral irradiance difference alone is 4.6%. Other factors such as the Libya-4 surface reflectance, 
gaseous absorption difference, and the absolute calibration of the two sensors will also impact the 
ratio. Bhatt et al. [29] shows that after accounting for the spectral differences using the SCanning 
Imaging Absorption SpectroMeter for Atmospheric CHartographY (SCIAMACHY) hyper-spectral 
radiances observed over Libya-4, the VIIRS M5 to MODIS Band 1 radiance ratio is nearly 1.0, which 
suggests that the absolute radiometric calibration of MODIS and VIIRS are consistent. The full 
investigation of the absolute radiance difference between other VIIRS and similar MODIS bands, 
mainly due to spectral disparity, is beyond the scope of this study. For bands where the MODIS and 
VIIRS spectral differences are much larger, the band dependent atmospheric absorption effects and 
BRDF differences may increase the temporal noise and possibly impact the trend analysis, which is not 
investigated in this study. The authors, however, believe that because the maintained sun-synchronous 
satellite orbital ground tracks of VIIRS repeat annually, any uncorrected bidirectional effects will also 
be replicated annually, and will not artificially alter the VIIRS stability analysis. A trend of the VIIRS 
ratios is then established and normalized to 1.0 at the midpoint of the VIIRS record. 

Figure 2. (a) Linear BRDF regression models derived from near-nadir TOA radiance 
measurements from Aqua-MODIS band 1 over Libya-4 (2002-2013) for forward (in red) 
and back scattering directions (in green); and (b) the VIIRS M5 observed divided by the 
predicted BRDF modeled in (a) radiance ratios before the normalization of the ratios over 
the VIIRS record. 
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3.2. DCC Approach 

Deep convective clouds are invariant targets located over the tropical tropopause layer (TTL) and 
are the brightest of all terrestrial invariant targets [15,30,31]. They have a high signal-to-noise ratio 
over visible wavelengths, exhibit nearly flat visible reflection spectra, and minimally impacted by 
water vapor and aerosols in TOA measurements [30,31]. Since DCC are the coldest targets at the 
tropics, they can be easily identified using a simple IR temperature threshold. Another noteworthy 
feature of DCC as invariant target is that they are viewed by (and thus applicable to) all satellite 
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sensors. The first attempt of using the DCC pixel-level radianees for ealibration purpose was made by 
Hu et al. [14] to examine the radiometrie stability of the CERES shortwave measurements on board the 
Tropical Rainfall Measuring Mission (TRMM), as well as of the Visible and Infrared Scanner (VIRS) 
and MODIS imagers on TRMM and Terra satellites, respectively. A similar approach was followed by 
Doelling et al. [15] to monitor the calibration drift of the Advanced Very High Resolution Radiometer 
(AVHRR) sensors on board the NOAA-16 and NOAA-17 satellites. Eater, Doelling et al. [30-32] 
made significant contributions towards refining the technique by analyzing the natural variability of 
the DCC response geographically, temporally, seasonally, and for differences of convection initiated 
over land and ocean, using 10-years of Aqua-MODIS data. The spectral characteristics of DCC were 
also examined using SCIAMACHY -measured hyper-spectral radiances and were reported to be flat in 
the VIS and NIR region [31]. 

In this study, the VIIRS DCC pixels are acquired over a fixed Tropical West Pacific (TWP) domain 
extending from 120°E to 160°E longitude, and 20°N to 20°S latitude, which has the highest frequency 
of DCC events. DCC radiance pixels within the TWP DCC domain are identified using an 1 1 pm 
brightness temperature (BT) less than 205.0° K threshold. To take advantage of the most isotropic part 
of the DCC reflectivity the identified DCC pixels are then filtered by SZA < 40°, VZA < 40°, 
10° < RAA< 170° [31]. To remove optically thin clouds, the following spatial standard deviation 
filters are used; obt< T0° K, and ovis< 3%, where OBiand avis are the standard deviations of the 
surrounding 3x3 pixel block in IR (11 pm) and VIS (0.65 pm) channels, respectively. The fdtered 
DCC radiance pixels are anisotropically corrected using the angular distribution model (ADM) from 
Hu et al. [14], which normalizes all DCC response to a common set of angular conditions. The normalized 
DCC pixel radiances are then compiled into monthly probability distribution functions (PDFs). 

Figure 3. (a) The VIIRS M5 channel DCC monthly PDFs over the TWP domain; 

and (b) the corresponding PDF mode and mean radiance trends showing the stability of the 

sensor. The ADM normalized radiance units are in W-m^^-sr ’-pm *. 
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Figure 3a shows the monthly PDFs of the normalized DCC radianees from the VIIRS M5 ehannel. 
The peaks of the PDFs are eonsistent as expeeted from a stable instrument. Over 200,000 DCC pixels 
were identified every month using this teehnique. The statistieal mean and mode are eomputed for eaeh 
monthly PDF, and are traeked over time to monitor the stability of the instrument. Figure 3b shows the 
linear regression of the monthly PDF modes and means to derive the stability of the VIIRS M5 ehannel, 
whieh is found to be 0.39%/yr using the mean and -0.16%/yr using the mode. It is noted that the mode 
statistie exhibits a third of the month-to-month variability than that of the mean statistie for this band. For 
VIS bands, the mode statistie is used for stability analysis. For SWIR bands, the mean statistie is used, 
sinee it provides more eonsistent monthly TO A radianees than the mode statistie [31]. 

4. Target Temporal Variability 

The monthly averaged observed desert ratios and monthly DCC anisotropie normalized radianees 
are further normalized to 1 at the midpoint of the VIIRS reeord and simple temporal linear regressions 
are applied to them. The results are shown in Figure 4 along with assoeiated statisties for quantitative 
analysis. The DCC natural variability is half of the desert natural variability for most VIS bands. For 
the 1.6 pm bands, the desert response is more stable, and for the Mil and B7 bands, the desert and 
DCC response stability is eomparable. The seasonal variations in the desert TOA observations are 
evident. At first glanee, no trends exeeeded 1%/yr and there were no band outlier trends. The 
desert-based trends were positive whereas the DCC derived trends were negative. The VIIRS 
ealibration drift is smaller, given that the drift is embedded in the invariant target natural variability. 
For deserts, the natural variability is the eombination of seasonal atmospherie absorption and aerosol 
effeets [20], and any short-term (e.g., owing to rain events) surfaee refleetivity ehanges [33]. For DCC, 
the natural variability is manifested in the temporal and spatial eoverage fluetuations during the annual 
migration of DCC with the sun [31]. This is espeeially true for the SWIR bands, where the signal to 
noise ratio is mueh smaller than for the VIS bands, and the refleetanee is dependent on the eloud 
mierophysies, whieh ehanges during the DCC eonveetive eyele. However, the NPP sun-synehronous 
orbit guarantees eonsistent loeal time sampling of the DCC eonveetive eyele. 

Beeause the VIIRS data reeord may be short for reliable trend analysis, the stability assessment is 
performed by eomparing the target variability against MODIS observations. If the temporal standard 
deviation of the VIIRS observations exeeeds those of the MODIS observations for similar bands 
during the VIIRS reeord, then the inereased observational noise ean be attributed to the VIIRS 
instrument eharaeteristies, sueh as ealibration and speetral band differenees. Then it follows that if the 
standard deviations of the VIIRS and MODIS observations are eonsistent during the VIIRS reeord, 
then both instruments have similar ealibration stability. It may also be the ease that the invariant target 
variability greatly exeeeds the ealibration uneertainty of both sensors. In this ease, the observed 
variability in the VIIRS ehannel response is merely due to the stability limits of the targets used. In 
order to deteet any short-term invariant target TOA radianee anomalies, the MODIS observations 
during the VIIRS time period ean be eompared to the full MODIS reeord. Not aeeounting for 
short-term invariant target brightening or darkening events may be interpreted as ealibration issues. 


Remote Sens. 2014, 6 


2817 


Figure 4. The normalized desert and DCC responses time-series for VIIRS RSB during the 
2012-2013 time period. 
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The same desert and DCC methodologies deseribed in this paper are performed on the MODIS 
Collection 6 radiances for stability analysis and the results are summarized in Table 2. Table 2 displays 
the temporal standard deviation of Aqua-MODIS invariant target radiances over the full record (July 
2002 to October 2013) and the VIIRS time frame (Febmary 2013 to October 2013), and the 
corresponding NPP VIIRS radiance stability. The Terra-MODIS record (January 2002 to October 
2013) is used to compare the 1.6 pm bands, due to inoperable detector issue with the corresponding 
Aqua-MODIS band [28]. All VIIRS and corresponding MODIS band standard deviations during the 
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VIIRS time frame are within 0.2%, except for the Mil SWIR band, where the VllRS standard 
deviation is 1.6% for both desert and DCC invariant targets. This is a reduction of 0.5% and 1.5% for 
deserts and DCC, respectively, for the Mil band, and may simply be due to the non-overlapping 
spectral bands. The VllRS M5 and 11 (0.65 pm) bands show a consistent decrease for both desert and 
DCC targets between 0.1% and 0.2% over Aqua-MODIS band 1. In fact the lowest standard deviation 
recorded is the VllRS M5 band using DCC at 0.4%. The VllRS 13 and MIO (1.6 pm) band shows an 
increase of 0.4% and 0.2% over Terra-MODIS over desert. This result is not confirmed by DCC, 
because the DCC SWIR band natural variability is much higher. To conclude, the VllRS response over 
Libya-4 and DCC in bands M3 (0.48 pm), M4 (0.55 pm), and M7 (0.86 pm) is stable on order of 
corresponding MODIS band, whereas the M5 (0.65 pm), 11 (0.65 pm), and Mil (2.25 pm) band 
responses are more stable than MODIS, and the 13 and MIO (1.61 pm) band responses are less stable 
than MODIS. The measured disparities between the VllRS and MODIS temporal standard deviations 
are probably due to differences in the instrument characteristics, such as SRF, that affect measured 
radiances over these targets. 


Table 2. Comparison of the observed desert (Libya-4) and DCC (TWP domain) target 
stability (%) computed from the full MODIS record (July 2002 to October 2013) and during 
the VllRS time frame (February 2013 to October 2013), and from VllRS observations. 


VIIRS/MODIS 

Band 


Libya-4 



DCC 


Observed MODIS Stability 
(%) 

Observed 

VIIRS 

Stability 

(%) 

Observed MODIS Stability 
(%) 

Observed 
VIIRS 
Stability (%) 

Full Reeord 

VIIRS 

Timeframe 

Full Reeord 

VIIRS 

Timeframe 

M3/B3 

1.1 

1.3 

1.4 

0.8 

0.5 

0.5 

M4/B4 

1.0 

1.2 

1.1 

0.7 

0.5 

0.5 

M5/B1 

0.9 

1.0 

0.8 

0.9 

0.6 

0.4 

M7/B2 

1.0 

1.0 

0.9 

MODIS band saturates 

0.5 

M10/B6 

0.8 

0.9 

1.1 

2.1 

2.1 

2.2 

M11/B7 

1.9 

2.1 

1.6 

3.9 

3.1 

1.6 

11/Bl 

0.9 

1.0 

0.8 

0.9 

0.6 

0.5 

I3/B6 

0.8 

0.9 

1.3 

2.1 

2.1 

2.2 


Figures 5 and 6 show the desert and DCC, respectively, TOA normalized radiances over the entire 
MODIS record. The VllRS time period is to the right of the red vertical dotted line. The standard 
deviations of the entire MODIS and VllRS time period are posted in both Figures 5 and 6 as well as in 
Table 2. Any short-term fluctuations over the invariant targets during the VllRS time period can be 
identified by visual examination. The MODIS desert standard deviations were greater by 0.1% to 0.2% 
during the VllRS time period than during the entire MODIS record for all VIS and SWIR bands, 
except for the M7 (0.86 pm) band, which was similar. This is easily identifiable in Figure 5, where the 
MODIS response slightly went down in most bands during 2011-2012, and recovered in early 2013. 
This could be the cause for the positive observed VllRS trends over Libya-4. 

On the contrary, the measured MODIS DCC standard deviations in VIS bands were smaller by 
0.2% to 0.3% during the VllRS time period, and are comparable to the corresponding VllRS bands 
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DCC standard deviations. The exception is the MODIS Band 7, where the DCC response was found to 
be less stable (3.1%) than the VIIRS Mil band (1.6%), and may be owed to the different spectral 
characteristics of these two bands. In addition, the detector striping in VIIRS-TEB is significantly 
better than the corresponding MODIS bands [8]. Improved RSB and TEB pixel co-registration and 
reduced noise improves the results of the DCC methodology. 


Figure 5. The normalized desert responses for MODIS RSB during the 2002-2013 time 
period. The standard deviation (stability) for the whole record and the VIIRS time frame is 
labeled in black and red (*) respectively. The dotted horizontal red line demarcates the 
beginning of the VIIRS record. 
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Figure 6. Same as Figure 5, except for the DCC invariant target. The MODIS band 2 
saturates for DCC target. 
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The DCC natural variability is half of the desert natural variability for most VIIRS VIS bands, for 
the 1.6 pm bands, the desert is more stable, and for the Mil band the stability is comparable. This is 
because the DCC spectrum is flat and has a nearly Lambertian reflectance over the VIS wavelengths, 
whereas the desert is impacted by seasonal impact of water vapor and aerosol variations. For the SWIR 
bands, the desert is more robust, because the DCC pixels are less reflective in this spectral domain and 
are dependent on the cloud microphysical properties, such as particle size, which undergo changes 
during the convection cycle. 
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5. Trend Significance 


If the invariant target observations represented only the VIIRS calibration drifts in Figure 4, then 
simple statistical trend significance tests can be performed. However, in this case the calibration trend 
is only significant if it exceeds the natural variability of the invariant target. The VIIRS linear 
regressions are evaluated using the statistical approach of trend detection as described by 
Weatherhead et al. [34,35]. This approach relates the minimum number of years of data required to 
confidently detect a given specified trend within the magnitude of the natural variability of a 
time-series of data. Or it can also be used to determine the minimum trend that can be detected with 
confidence based on the record length. 

According to Weatherhead et al [35], the following equation can be used to compute the required 
number of years (N) to detect a trend of magnitude 0o at the 95% confidence level and with a 
50% probability. 


N = (- 


l^^ol 


jl + 0 

Jl-0 


y/3 


( 1 ) 


where, ajv and 0 are the month-to-month variability and 1 -month lag autocorrelation in the data 
time-series, respectively. If the observations are auto-correlated or not independent from each other, 
then trend detection will take longer. In this study, the equation is slightly modified to use N as an 
input and solve for coo, which now gives the minimum detectable trend for a given data time-series. It 
should be noted that the units of on and coo are the same in this equation, on is expressed in percentage 
by dividing it by the mean value. Similarly, N is computed in terms of years by dividing the total 
number of months by 12. 


Table 3. Comparison of the observed VIIRS desert and DCC trends (%/yr) with the 
minimum detectable trend computed at 95% confidence level and with 50% probability 
using Equation (1). 


VIIRS 

Band 


Libya-4 


DCC 

Observed Trend 
(%/yr) 

Minimum Deteetable 
Trend ± (%/yr) 

Observed Trend 
(%/yr) 

Minimum Deteetable 
Trend ± (%/yr) 

M3 

0.77 

1.83 

-0.08 

0.53 

M4 

0.50 

1.66 

-0.67 

0.82 

M5 

0.71 

1.34 

-0.16 

0.54 

M7 

0.63 

1.70 

-0.11 

0.44 

MIO 

0.71 

2.00 

-0.87 

3.24 

Mil 

0.05 

2.60 

-0.51 

2.2 

11 

0.31 

1.26 

-0.42 

0.52 

13 

0.76 

2.50 

-0.87 

3.22 


Table 3 displays the band-specific observed band trend shown in Figure 4 and the associated 
minimum detectable trend (giq), based on Equation (1), which can be confidently detected given the 
natural variability of the data time series. All trends have been expressed in %/yr to facilitate the 
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comparison of the two trends. The comparison shows that all of the observed trends are within the 
eorresponding minimum deteetable trend amounts, thus indieating that the observed trends do not rise 
above the natural variability. It ean also be stated that the VIIRS ealibration drift does not exeeed 
minimum deteetable trend. The VIIRS VIS band ealibration drift is within 0.54%/yr, exeept for the M4 
band (0.55 pm) based on the DCC invariant target. The VIIRS band ealibration drift is contained 
within 2.5%/yr. 

Generally, the magnitude of the natural variability ean be estimated by the standard deviation of the 
monthly TOA observations. In addition, the standard deviation is eonstant in time (Figures 5 and 6) 
after the reeord reaehes a few years. In this ease, the VIIRS reeord eontains a short-term desert 
fluetuation, and the early MODIS DCC record contains noisy pixel-level radianees. If one assumes that 
the future Libya-4 VIIRS lifetime reeord will resemble the MODIS full reeord, and that the future 
DCC VIIRS lifetime reeord will resemble the current VIIRS time frame, one can compute the number 
of years required to deteet a trend with a magnitude equal to the temporal standard deviation 
per deeade. Beeause the natural variability or monthly standard deviation is band-dependent, the 
number of years needed to deteet a trend based on the overall standard deviation should be similar 
among bands. 

Table 4 displays the number of years required to significantly detect a trend with a magnitude equal 
to the standard deviation per deeade. For Libya-4, the number of years required is ~11 years, and for 
DCC ~9 years. If the DCC VIS band trends are doubled the number of years is redueed to ~5 years 
(not shown). This supports that any VIIRS ealibration drifts greater than 1%/deeade can be deteeted in 
~5 years using DCC targets and in ~1 1 years using Libya-4. 


Table 4. The band-speeifie trends are arbitrarily based on the invariant target standard 
deviations divided by a deeade. The number of years needed to deteet the given trend at the 
95% confidenee level and with 50% probability using Equation (1), where the Libya-4 time 
series is based on the MODIS full reeord, and the DCC is based on the VIIRS reeord is 
also listed. 


Libya-4 DCC 

(MODIS Full Record) (VIIRS) 

Band Trend ± (%/Decade) Years To Detect Trend Trend ± (%/Decade) Years to Detect Trend 


M3 

1.1 

11 

0.5 

8 

M4 

1.0 

11 

0.5 

10 

M5 

0.9 

10 

0.4 

9 

M7 

1.0 

11 

0.5 

7 

MIO 

0.8 

14 

2.2 

9 

Mil 

1.9 

10 

1.6 

9 

11 

0.9 

10 

0.5 

7 

13 

0.8 

17 

2.2 

9 


6. Concluding Remarks 

This study has provided an initial assessment of the radiometrie ealibration stability of the S-NPP 
VIIRS RSB of the NASA Land PEATE AS3100 produet arehived at ASDC using the Eibya-4 desert 
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and the DCC pixels loeated over the TWP region as two independent radiometrieally stable targets. 
This VIIRS produet has removed all known ealibration anomalies, sueh as the sensor degradation of 
the UV bands and the bands impaeted by tungsten eontamination [8]. The stability of the VIIRS 
multispeetral imager is important for the CERES projeet to derive temporally eonsistent and reliable 
eloud retrievals. The VIIRS stability analysis was limited by the natural variability of the invariant 
targets and the brief VIIRS observational reeord. For VIS bands, the DCC target has half of the 
temporal variability eompared to Eibya-4, whereas for the SWIR bands, Eibya-4 has less temporal 
variability. The desert TOA BRDF normalized radianee reveals a seasonal dependeney due to the 
water vapor and aerosol fluetuations. The low DCC refleetanee and eloud mierophysieal variations 
dependent on the eonveetive eyele inerease the temporal variability in the SWIR bands. The temporal 
standard deviation of monthly VIIRS DCC TOA observations is within 0.53% for VIS bands and for 
SWIR bands the temporal standard deviation of VIIRS Eibya-4 TOA observations is within 1.56%. 
The Eibya-4 target stability is hampered by a brightening event that oeeurred during the VIIRS reeord. 
The improved VIIRS DCC VIS band temporal stability eompared with MODIS eonfirms that the 
VIIRS deteetor noise, or striping, is signifieantly smaller than that of MODIS [8]. Otherwise the VIIRS 
temporal stability is similar to the eorresponding MODIS band temporal stability. 

The VIIRS ealibration drift is embedded in the invariant target natural variability. Only drifts that 
exeeed the natural variability ean be deteeted using these targets. The minimum deteetable trend was 
eomputed over the VIIRS observed reeord, along with linear regressions based on the measurements. 
The desert measured trends ranged from 0.05%/yr to 0.77%/yr. The DCC trends ranged from 
-0.08%/yr to -0.87%/yr. None of these trends exeeeded the minimum deteetable trend based on the 
observed natural variability over the invariant targets. Only longer time reeords ean eonfirm smaller 
VIIRS ealibration trends using these targets. The faet that the desert and DCC trends are of opposite 
sign reveals that the target variability masks the true VIIRS ealibration signal. For VIS bands, it is 
predieted that ~9 years are required to deteet a VIIRS ealibration trend greater than 0.55%/deeade 
using the DCC target and ~5 years to deteet a trend of 1%/deeade. For SWIR bands ~II years are 
required to eonfirm a trend that exeeeds ~ 1%/deeade. 

Currently, the VIIRS instrument ealibration stability based on the two stable Earth targets diseussed 
in this paper is found similar to or better than the MODIS instrument. This study also supports that the 
NPP VIIRS AS3I00 subsetted data from NASA Eand PEATE should provide eonsistent eloud 
retrievals for the CERES FM-5 flux data produets during the time of study. It must be noted that any 
sean mirror ealibration drifts for sean positions outside of 40° (DCC) and 10° (Eibya-4) VZA were not 
evaluated in this study. 
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